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Abstract 



The 7]ir~ system has been studied in the reaction ir~p — > r]ir~p at 18 GeV/c. 
A large asymmetry in the angular distribution is observed indicating interfer- 
ence between L-even and L-odd partial waves. The (1320) is observed in the 
J PC = 2++ wave, as is a broad enhancement between 1.2 and 1.6 GeV/c in 
the 1 _+ wave. The observed phase difference between these waves shows that 
there is phase motion in addition to that due to ai (1320) decay. The data can 
be fitted by interference between the 02(1320) and an exotic 1 _+ resonance 

with M = (1370 ±16 +3°) MeV/c 2 and T = (385 ±40 MeV/c 2 . 
13.60.Le, 13.85.Fb, 14.40.Cs 
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The question of whether or not hadrons outside the scope of the constituent quark 
model exist is one whose answer speaks directly to the fullness of our understanding of 
quantum chromodynamics (QCD) |L|. However, non-gg mesons (or exotic mesons) have 
proven difficult to distinguish from the many conventional qq states which populate the 
various mesonic spectra. For this reason, much attention has been focused on those states 
with manifestly exotic J PC quantum numbers. 

A qq meson with orbital angular momentum £ and total spin s must have P = (— 
and C = (— l) e+s . Thus a resonance with J PC = , + ", 2 + ~, ... must be exotic. 
Such a state could be a gluonic excitation such as a hybrid (qqg) or glueball (2g, 3g, ...), or 
a multiquark (qqqq) state. In a relative P wave (L=l), the r]ir~ system has J PC = 
Having isospin 1=1, it could not be a glueball, but it could be a hybrid or a multiquark 
state. 

Production and decay properties of exotic states have been predicted using several models 
0-^. A calculation based upon the MIT bag model predicts |J that a 1~ + hybrid (qqg) 
will have a mass near 1.4 GeV/c 2 . On the other hand, the flux-tube model predicts the 
mass of the lowest-lying hybrid state to be around 1.8 GeV/c 2 . Characteristics of bag-model 
S-wave multiquark states (which would have J p = + , 1 + , or 2 + ) have been predicted ]7[ 
but those for a 1~ state have not. Finally, recent lattice calculations |§ of the 1~ + hybrid 
meson estimate its mass to be in the range of 1.7 to 2.1 GeV. 

The 7]ir system has been studied in several recent experiments, with apparently inconsis- 
tent results. Aide et al. 0, in a study of 7r~p interactions at 100 GeV/c at CERN (the GAMS 
experiment), claimed to observe a 1 _+ state in the rjn system at 1.4 GeV/c 2 produced via 
unnatural parity exchange (the Po partial wave — the naming convention is discussed be- 
low) |H| . Aoyagi et al. |TTJ, in a ir~p experiment at 6.3 GeV/c at KEK, observed a rather 



narrow enhancement in the 7]n~ system at 1.3 GeV/c 2 in the natural parity exchange 1~ + 
spectrum (P+). Beladidze et al. ||12|| , in the VES experiment at IHEP, (tt~N interactions 
at 37 GeV/c) also reported a P + signal in the t]ir~ state, but their signal was broader and 
had a significantly different phase variation from that of the KEK experiment. While the 



phase difference between the P + and D + waves was independent of rjn mass in the KEK 
analysis, that phase difference did show significant mass dependence in the VES analysis. 
(Since the phase variation for the D + wave follows a classic Breit-Wigner pattern for the 
02(1320) meson, the phase difference between these waves can determine the phase variation 
of the unknown P + wave.) 

Here we study the rjir^ system in the reaction ir~p — > r]7T~p at 18 GeV/c. Our data 
sample was collected in the first data run of E852 at the AGS at Brookhaven National 
Laboratory with the Multi-Particle Spectrometer (MPS) |L3f using a liquid hydrogen target. 
The MPS, which was equipped with six drift-chamber modules [14| and three proportional 



wire chambers, was augmented by: a four-layer cylindrical drift chamber surrounding the 



target [|T^]; a soft-photon detector consisting of 198 blocks of thallium-doped cesium iodide 
fL6|| also surrounding the target; a window- frame lead-scintillator photon- veto counter; a 
large drift chamber; and a 3045-element lead-glass detector (LGD) [17[ downstream of the 



MPS. Further details are given elsewhere [18| . 

A total of 47 million triggers which required one forward-going charged track, one recoil 
charged track, and an LGD trigger-processor signal enhancing high electromagnetic effective 
mass was recorded. Of these, 47,200 events were reconstructed which were consistent with 
the rjii~p (77 — > 27) final state. These events satisfied topological and fiducial volume cuts, 
as well as energy/momentum conservation for production and for the 77 — ^ 2^ decay with a 
confidence level > 10% ||19|| . The 27 mass resolution at the 77 mass is a = 0.03 GeV/c 2 . 



The 02(1320) is the dominant feature of the t]ti~ mass spectrum shown in Fig. |T|a. Back- 
ground has been estimated using side bands in both the 2-7 mass distribution and the 
missing-mass distribution, thus taking into account background from non-77 sources as well 
as from sources due to production of other final states. The background level is approxi- 
mately 7% at 1.2 GeV/c 2 , falling to 1% at 1.3 GeV/c 2 . 

The acceptance-corrected distribution of \t'\ = \t\ — |t| m i n , where t is the the four- 
momentum-transfer, is shown for \t'\ > 0.08(GeV/c) 2 in Fig. |IJb. (Our acceptance is quite 
low below 0.08 (GeV/c) 2 due to a trigger requirement.) The shape of this distribution is con- 



sistent with previous experiments and has been shown to be consistent with natural-parity 



exchange production in Regge-pole phenomenology |2"0| , |21" |. 

The acceptance-corrected distribution of cos#, the cosine of the angle between the rj and 
the beam track in the Gottfried- Jackson frame |2!| of the rjTr" system, is shown in Fig. 0a 
for 1.22 < M(i]tt~) < 1.42 GeV/c 2 . There is a forward-backward asymmetry in cos#. The 
asymmetry for | cos#| < 0.8 is plotted as a function of r]7T~ mass in Fig.^b. The asymmetry 
is large, statistically significant and mass dependent. Since the presence of only even values 
of L would yield a symmetric distribution in cos 8, the observed asymmetry requires that 
odd-L partial waves be present to describe the data. 



A partial-wave analysis (PWA) |23|j24| based on the extended maximum likelihood 



method has been used to study the spin-parity structure of the r]ir~ system. The partial 
waves are parameterized in terms of the quantum numbers J PC as well as m, the absolute 
value of the angular momentum projection, and the reflectivity e (which is positive (neg- 
ative) for natural (unnatural) parity exchange ||25|| ). In our naming convention, a letter 
indicates the angular momentum of the partial wave in standard spectroscopic notation, 
while a subscript of means m = 0, e = — 1, and a subscript of +(— ) means m = 1, e 
= +1(— 1). Thus, So denotes the partial wave having J PC m e = ++ 0~, while P_ signifies 
1~ + 1~, D + means 2 ++ l + , and so on. We consider partial waves with m < 1, and we assume 
that the production spin-density matrix has rank one. 

The experimental acceptance is determined by a Monte Carlo method. Peripherally- 



produced events are generated |26| with isotropic angular distributions in the Gottfried- 



Jackson frame. After adding detector simulation p7| , the Monte Carlo event sample is 
subjected to the same event-selection cuts and run through the same analysis as the data. 
The experimental acceptance is then incorporated into the PWA by using these events to 
calculate normalization integrals (see ref . |23| ) . 

Goodness-of-fit is determined by calculation of a x 2 from comparison of the experimental 
moments with those predicted by the results of the PWA fit. A systematic study has been 
performed to determine the effect on goodness-of-fit of adding and subtracting partial waves 



of J < 2 and m < 1. All such waves have been included in the final fit. We have also 
performed fits including partial waves with J = 3 and J = 4. Contributions from these 
partial waves are found to be insignificant for M(r/7r _ ) < 1.8 GeV/c 2 . Thus, PWA fits 
shown or referred to in this letter include all partial waves with J < 2 and m < 1 (i.e. 
S , Po, D , D_, P + , and D + ). The background described above was included as a 
non-interfering, isotropic term of fixed magnitude. 

The results of the PWA fit in 40 MeV/c 2 bins for 0.98 < M^tt") < 1.82 GeV/c 2 and 
0.10 < \t\ < 0.95 GeV 2 are shown in Fig.[|a-c. Here, the acceptance-corrected numbers 
of events predicted by the PWA fit for the D + and P + waves and their phase difference 
A$(D + — P + ) are shown as a function of M(r]7r _ ). There are eight ambiguous solutions in 
the fit p^f28p9"f , each of which leads to the same angular distribution. We show the range of 



fitted values for these ambiguous solutions in the vertical rectangular bar at each mass bin, 
and the maximum extent of their errors is shown as the error bar. The 0,2(1320) is clearly 
observed in the D + partial wave (Fig.|3]a). A broad peak is seen in the P + wave at about 1.4 
GeV/c 2 (Fig. 0b). A$(D + — P + ) increases through the 02 (1320) region, and then decreases 
above about 1.5 GeV/c 2 (Fig.|3]c). The intensities for the waves of negative reflectivity (not 
shown) are generally small and are all consistent with zero above about 1.3 GeV/c 2 . 

These results are quite consistent with the VES results [|l^]. In particular, the shape of 
the phase difference is virtually identical to that reported by VES. (The magnitude of the 
phase difference is shifted by about 20° relative to that of VES.) 

Consistency checks and tests of the data have been carried out to determine whether the 
observation of the structure in the P + wave could be an artifact due to assumptions made in 
the analysis or to acceptance problems. These include: fitting the data in restricted ranges of 
the decay angle; inclusion of higher angular momentum states; fitting the data with various 
t cuts; fitting the data using different parametrizations of the background; making cuts on 
other kinematic variables such as the ir~p or the rjp effective masses; and fitting data using 
events with r] — > 7r + 7r~7r° decays (with rather different acceptance from the 27 events). The 
results are very stable and, in particular, the behavior of A$(D + — P + ) does not change in 
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any of these checks. 

Fits were also carried out on Monte Carlo events generated with a pure D + wave to 
determine whether P + -wave structure could be artificially induced by acceptance effects, 
resolution, or statistical fluctuations. We do find that some P + intensity can be induced by 
resolution and/or acceptance effects. Such "leakage" leads to a P + wave that mimics the 
generated D + intensity (and in our case would therefore have the shape of the ei2(1320)) 
with a A$(D + — P + ) that is independent of mass. Neither property is seen in our result. 

In an attempt to understand the nature of the P+ wave observed in our experiment, 
we have carried out a mass- dependent fit to the results of the mass-independent amplitude 
analysis. The fit has been carried out in the rjn mass range from 1.1 to 1.6 GeV/c 2 . The 
input quantities to the fit included, in each mass bin, the P + -wave intensity; the D + -wave 
intensity; and the D + — P + phase difference. Each of these quantities was taken with its 
error and correlation coefficients from the result of the amplitude analysis. In this fit, we 
have assumed that the D + -wave and the P + -wave decay amplitudes are resonant and have 
used relativistic Breit-Wigner forms JJIJ for these amplitudes. We introduce a constant 
relative production phase between the P + -wave and D + -wave amplitudes. The parameters 
of the fit included the D + -wave mass, width and intensity; the P + -wave mass, width and 
intensity; and the D + — P + production phase difference. One can view this fit as a test of 
the hypothesis that the correlation between the fitted P-wave intensity and its phase (as a 
function of mass) can be fit with a resonant Breit-Wigner amplitude. 

Results of the fit are shown as the smooth curves in Fig. ^|a, b, and c. The mass and 
width of the J PC = 2++ state (Fig. 0a) are (1317 ±1 ±2) MeV/c 2 and (127 ±2 ±2) MeV/c 2 
respectively fl3~I] . (The first error given is statistical and the second is systematic f32|.) The 
mass and width of the J PC = 1~ + state as shown in Fig.^b are (1370 ±16 ^) MeV/c 2 
and (385 ±40 ^) MeV/c 2 respectively. Shown in Fig.|d are the Breit-Wigner phase 
dependences for the a 2 (1320) (line 1) and the P + waves (line 2); the fitted D + — P + production 
phase difference (line 3); and the fitted D + — P + phase difference (line 4). (Line 4, which is 
identical to the fitted curve shown in Fig.||]c, is obtained as line 1 — line 2 + line 3.) 
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The fit to the resonance hypothesis has a % 2 /dof of 1.49. The fact that the production 
phase difference can be fit by a mass-independent constant (of 0.6 rad) is consistent with 
Regge-pole phenomenology in the absence of final-state interactions. If one fits the 
data to a non-resonant (constant phase) P + wave, and also assumes a Gaussian intensity 
distribution for the P + wave, one obtains a fit with a x 2 /dof of 1.55. In this case, the 
observed phase dependence on mass is attributed to a rapidly varying production phase 



3J|. Such a phase variation cannot be excluded, but is not expected for any known model. 



Note that for this non-resonant hypothesis one must have a separate hypothesis for the 
observed structure in the P + intensity — a structure which is explained naturally by the 
resonance hypothesis. We thus conclude that there is credible evidence for the production 
of a J PC = 1 h exotic meson. 
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Without their outstanding efforts, the results presented here could not have been obtained. 
We would also like to acknowledge the invaluable assistance of the staffs of the AGS and 
BNL, and of the various collaborating institutions. This research was supported in part 
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FIGURES 




FIG. 1. a.) The rjir effective mass distribution, b.) Distribution of \t'\ = \t\ — \t\ 
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FIG. 2. Distributions of a.) the cosine of the decay angle in the Gottfried- Jackson frame for 
events with 1.22 < M(r]ir~) <1.42 GeV/c 2 , and b.) the forward-backward decay asymmetry as a 
function of M{t]tt~). The asymmetry = (F — B)/(F + B) where F(B) is the number of events for 
which the 77's momentum is forward (backward) in the Gottfried-Jackson frame. The dashed curve 
and the right-hand scale in a.) show the acceptance in this mass region. 
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FIG. 3. Results of the partial wave amplitude analysis. Shown are a.) the fitted intensity- 
distributions for the D + and b.) the P+ partial waves, and c.) A<E>(D + — P + ), their phase difference. 
The range of values for the eight ambiguous solutions is shown by the central bar and the extent 
of the maximum error is shown by the error bars. Also shown as curves in a.), b.), and c.) are 
the results of the mass dependent analysis described in the text. The lines in d.) correspond to 
(1) the fitted D+ Breit-Wigner phase, (2) the fitted P+ Breit-Wigner phase, (3) the fitted D+-P+ 
relative production phase, and (4) the overall D + — P + phase difference as shown in c.) but with a 
different scale. 
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